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ABSTRACT 

After a long quiescence of three decades, the transient X-ray pulsar 4U 1901+03 
became highly active in 2003 February. From the analysis of a large number of Rossi 
X-ray Timing Explorer/ Proportional Counter Array (RXTE/PCA) observations of 
this source, we report here the detection of X-ray flares, a broadening of the pulse 
frequency feature and Quasi Periodic Oscillations (QPOs). The X-ray flares showed 
spectral changes, had a duration of 100 s - 300 s, and were more frequent and stronger 
during the peak of the outburst. In most of the observations during the outburst we 
also detected a broadening of the pulse frequency peak. We have also found intensity 
dependent changes in the pulse profile at very short timescales. This reveals a coupling 
between the periodic and the low frequency aperiodic variabilities. In addition, near the 
end of the outburst we have detected a strong QPO feature centered at ^0.135 Hz. The 
QPO feature is broad with a quality factor of 3.3 and with an rms value of 18.5±3.1%. 
Using the QPO frequency and the X-ray luminosity during the QPO detection period 
we estimated the magnetic field strength of the neutron star as 0.31 xlO 12 G which is 
consistent with the value inferred earlier under the assumption of spin equilibrium. 
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1 INTRODUCTION 

Most of the accretion powered X-ray pulsars belong to the 
class of High Mass X-ray Binaries. The companion star may 
be a Be star or an OB supergiant. Most of the Be star sources 
are transients which are normally detected during outbursts. 
There are two kinds of outbursts, normal outburst (Type-I) 
and giant outburst (Type-II). Be stars have circumstellar 
disks and when the neutron star passes through the disk, it 
produces Type-I outbursts. The Type-II outbursts are pro- 
duced during the episodes of large mass outflow from the 
Be star. Transient X-ray pulsars are excellent candidates to 
investigate properties of the accretion process that are in- 
tensity dependent. 

The recurrent High Mass X-ray binary pulsar 4U 
1901+03 was discovered by Uhuru and Vela 5B satellite 
in 1970-1971 wh ich was the first detection of an outburst 
from this source |Forman et al.lll976l ; IPriedhorskv fc Terrell 
Il984l h A second outburst in 2003 was reported by Rossi X- 
ray Timing Explorer- All Sky Monitor (RXTE-ASM). After 
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the detection by ASM, detailed long term observations of 
this source was carried out with the Proportional Counter 
Array (RXTE-PCA). 

RXTE observations revealed pulsations of 2.763 s and 
the pulse frequency variation suggested an orbital period of 
22.58 days and eccentricity of about 0.035. No optical or IR 
coun terpart has been found for this source (|Gallowav et al.l 
l2005h . 4U 1901+03 is included in the unusual class of X-ray 
binaries with wide orbits (P or b ^ 20 days) but low eccen- 
tricity. The source was also observed in the hard X-ray band 
with the International Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) betw e en 2003 Ma r ch 10 and 2003 April 13 
(|Molkov et al.ll2003l hl Chen et all (|2008l ) carried out a pulse 
profile analysis of this source and found that the pulse pro- 
file is luminosity dependent. A phase resolved spectroscopic 
study revealed a stron g pulse phase dependence of the spec- 
trum (|Lei et alj|2009h . 

In the following sections we present a detailed timing 
analysis of RXTE/PCA archival data of 4U 1901+03 and 
report the discovery of several interesting features including 
flares, a broadening of the spin frequency peak, and a QPO 
at ~0.135 Hz. 
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Figure 1. The top panel shows the one day averaged ASM light 
curve and the bottom panel shows the background subtracted 
PCU2 light curve of 4U 1901+03 with a bin size of 200 s. The 
arrow marks indicate the period of QPO detection. 



2 OBSERVATIONS AND DATA REDUCTION 

The HMXB pulsar 4U 1901+03 was observed with the 
RXTE in 2003 during an outburst after a long interval of 
quiescence. The RXTE/PCA observations started on 2003 
February 10 and lasted for around 5 months. Though the 
RXTE/PCA has a large field of view of 1 square degree, 
any contamination from nearby sources is negligible because 
there are no bright X-ray sources within 1° from this object. 
There were a total of 68 pointings during this time period 
and data is available for a total useful exposure of 398 ks. 
The PCA instrument of RXTE consists of five Proportional 
Counter Units (PCU) covering an energy range of 2-60 keV 
with an effective area of 6500 cm 2 . The ASM is sensitive 
to X-ray photons between 1.5-12 keV. Table [1] gives a log 
of the observations. For our timing analysis we used data 
from all the PCA observations and during most of the ob- 
servations PCU0 and PCU2 were operational. The data re- 
duction was performed using the software package FTOOLS 
version 6.5. The background count rate was simulated using 
the tool pcabackest and was subtracted from the total count 
rate. Appropriate background models are used for the back- 
ground simulation. Caldb version 1.0.1 and pcarsp version 
11.7 were used for the creation of the energy response matrix 
of the PCA detections. 



3 ANALYSIS AND RESULTS 

3.1 Light curve and Flares 

Many accretion powered pulsars show aperiodic variability, 
like flares, with timescales of a few seconds to a few hundred 
seconds. Here we are presenting results from the analysis of 
data from the RXTE/PCA long term observations of 4U 
1901+03, in which we have detected a large number of X- 
ray flares. 

One day averaged ASM light curve is shown in the top 
panel of Figure [T] for a 250 day period around the outburst. 
The outburst started around MJD 52660 and reached a max- 
imum co unt rate of 18 s -1 in t he ASM detectors (about 240 
mCrab) l|Gallowav et al.ll2005h and then it declined linearly 
reaching a count rate of about 1 counts/sec at MJD 52840. 
We extracted light curves from all the 68 PCA observations 
with a time resolution of 0.125 s. The entire background sub- 
tracted PCU2 light curve from all 68 observations is shown 
in the bottom panel of Figure [T] This 2-60 keV light curve is 
generated with a time resolution of 200 s. The arrow marks 
in Figure 1 indicate the periods of QPO detection described 
later. 

In most of the observations carried out during the out- 
burst, 4U 1901+03 exhibited X-ray flares. The flares oc- 
curred randomly with different recurrence time between a 
few hundred to thousand seconds and flare duration of 100 
to 300 seconds. Pulsations are also clearly seen during the 
flares. Some of the flares had sharp rise and slow decay with 
multiple peak structure. In Figure [2] we have shown three 
light curves selected on the basis of intensity of flare activity 
and a fourth without flares but strong short term variability. 
The light curves are plotted with a bin size equal to the spin 
period of ~2.763 s so that the variability associated with the 
pulsation is omitted. At the peak of the flare, the luminosity 
is higher by a factor of upto 2.5. The flares are found to be 
more frequent and stronger near the peak of the outburst 
indicating a relation with the mass accretion rate. 

We calculated the hardness ratios (soft hardness ratio: 
6-10 keV / 2-6 keV and hard hardness ratio: 10-30 keV / 6- 
10 keV) of the light curves using the event mode data (Good 
Xenon mode) from the PCA. The hardness ratios of the data 
corresponding to panel (b) of Figure [2] are shown in Figure 
along with the summed light curve in the bottom panel. 
During the flares, the soft hardness ratio is found to increase 
while the hard hardness ratio is seen to be decreasing. This 
implies that the flares are stronger in the medium energy 
band which we have also found in the analysis of the flare 
spectrum described later. A detailed hardness ratio analysis 
showed that all the flares have similar energy dependence. 

3.2 Power density spectrum 

Quasi Periodic Oscillations (QPOs) have been reported in 
many Accretion Powered Pulsars, both HMXBs and LMXBs 
(see James et al. 2010 for a complete list of the sources). In 
the accreting high magnetic field X-ray pulsars, the QPO 
appears as a broad peak in the power density spectrum. 
We have carried out a timing analysis using data from all 
the RXTE-PCA observations to search for QPOs and other 
aperiodic features in this pulsar. 

Standardl mode data of RXTE/PCA with a time reso- 
lution 0.125 s was used for the timing analysis. We extracted 
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Figure 2. Some representative RXTE/PCA light curves of 4U 
1901+03. The binsize used here equals the spin period of the 
neutron star. The start time of the light curves are (a)52685 15:34 
(b)52685 21:57 (c)52689 13:08 (d)52686 18:25 MJD respectively 



> 

CD 



o 
I 

CD 
\ 
O 
CO 

I 



0.95 



0.9 



0.75 



„ 4000 



a 3000 

a 
a 

8 2000 



- 1 

w 


1 ' 


4 


-i — i — i — i — i— 


1 >• 

m 






^^^^ 




H 











1000 2000 
Time (s) 



3000 



Figure 3. Background subtracted RXTE/PCA light curves (Fig- 
ure 2b: 52685 21:57 MJD) and hardness ratios in two different 
energy band pairs are shown here. 



light curves from all the archival data and created Power 
Density Spectrum (PDS) using the FTOOL-powspec. The 
light curve was divided into stretches of length 1024 s and 
the power spectra obtained from 2-5 consecutive segments 
were averaged to produce a set of averaged PDS from which 
the expected white noise level was subtracted. The PDS were 
normalised such that normalised integral gives the squared 
rms fractional variability. 



For most of the observation period, the PDS was a fea- 
tureless continuum except for the pulse peak and its har- 
monics. The PDS reveals no QPO feature during the time 
of the outburst but in many of the observations we detected 
a broadening of the spin frequency peak. Near the end of the 
outburst, the broadening of the spin frequency peak disap- 
peared and a broad QPO feature at ~0.135 Hz appeared in 
the PDS. The three PDS with and without flares are shown 
in Figure U (top panel: PDS for the period of flare; middle 
panel: non flare; bottom panel: QPO detection). During the 
QPO detection we also detected a large pulsed fraction (upto 
60%) and also a large rms (18%) of the QPO. If these two 
features were to be from two different soruces in the FOV, 
then even if the pulsating source was to be 100% pulsed, the 
source with QPOs is required to have a rms of 45%. This is 
much larger than any QPO feature in any type of X-ray bi- 
nary. We, therefore, conclude that the QPO feature is indeed 
from the X-ray pulsar. The PDS with the QPO was fitted 
with a model consisting of a Power law and Lorentzian for 
the continuum and a second Lorentzian for the QPO feature. 
The frequency bins corresponding to the pulse peak and its 
harmonics seen in Figure [4] were ignored while fitting the 
continuum. Center frequency of the QPO feature was found 
to change from a value of 0.143 ± 0.003 Hz to 0.130 ± 0.006 
Hz in 3 days. The QPO feature has a width of about 0.04 
Hz and a quality factor of about 3.3. The significance of the 
detection of the QPO was different on different days varying 
from 6 a to 7 a. The background corrected rms value of the 
QPO feature is 18 ± 3%. We calculated the energy depen- 
dence of the QPO rms using the event mode data (Good 
Xenon mode) of the PCA and found that the QPO feature 
is not detectable above 11 keV. 

A broadening of the X-ray pul se frequency peak has 
been seen in several other sources jLazzati fc Stellal Il997l : 
iReig et al.ll2006l ; iRaichur fc Paulll2008l : I Jain et al.ll201Clh and 
is thought to be due to a coupling of the low frequency vari- 
ations with the pulsation. To investigate this in detail, we 
have created intensity dependent pulse profiles. Light curve 
from RXTE-PCA observations were first created with bin 
size same as the pulse period. Using the X-ray intensity of 
each single pulse, separate intensity dependent time win- 
dows were created and these windows were then used to ob- 
tain intensity dependent pulse profiles of 4U 1901+03. The 
intensity dependent pulse profiles obtained from two light 
curves, with and without flares (panel b & d in Figure 
are shown in the top two panels of Figure [5] It is to be noted 
that broadening of the pulse frequency peak is present in all 
the light curves shown in Figure [5] In bottom panel of Fig- 
ure [5] we have shown the intensity selected pulse profiles 
during a period when the light curve shows a QPO feature 
but no broadening of the X-ray pulse frequency peak. Note 
that the three PDS shown in Figure[4]and Figure [5]are from 
the same set of light curves (top: flare, middle: non-flare, 
and bottom: QPO). 

In the top two panels of Figure[5]one can clearly see that 
at low intensity, the pulse profile has two unequal peaks. At 
higher intensity, both the peaks become stronger, the dip be- 
fore the main peak disappears while the dip after the main 
peak becomes more clear. The intensity dependence of the 
pulse profile appears to be quite similar for the flaring and 
non-flaring episodes. From the bottom panel of Figure [5] we 
can see that the episode that does not show a pulse peak 
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Figure 4. Power density spectra with broadening (top and mid- 
dle panel) and without broadening of the pulse peak (bottom 
panel) are shown in this figure. 



broadening also does not have any significant intensity de- 
pendence of the pulse profile. 
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Figure 5. In each panel of this figure, pulse-profiles at three 
different intensity levels are shown. Top panel: when the light 
curve shows flares (Figure 2b); middle panel: when the light curve 
does not show flares (Figure 2d); bottom panel: when the light 
curves shows QPOs and no broadening of the pulse frequency 
peak (52825 18:09 MJD). 



Table 2. Best fit Spectral Parameters of 4U 1901+03 



3.3 Spectral Analysis 

We have done a spectral analysis of the source during the 
QPO detection, flaring and non flaring state to investigate if 
there are any differences. Standard2 data of the PCA were 
used for extraction of the energy spectrum in 2.5-25 keV 
energy band. All the three spectra were fitted with a model 
consisting of an absorption component, a power law with a 
high-energy cutoff and a gaussian emission line. A system- 
atic error of 0.5% was added to the first and second spec- 
tra to account for the calibration uncertainties while the 
non-flaring spectrum required a systematic error of 1%. The 
spectra, along with the respective best-fit models and the 
residuals are shown in Figure [5] (top panel: spectrum for the 
period of flare; middle panel: non flare; bottom panel: QPO 
detection). The best fitted spectral parameters are given in 
Table [2] along with the reduced %r values. The spectrum 
taken when QPOs are present and broadening of spin fre- 
quency peak is absent has a larger photon index and large 
value of the high energy cut-off compared to the spectrum 
when QPOs are absent and spin frequency peak shows a 
broadening, non flaring spectra gives Xr °f 1-3 and 1.6 re- 
spectively with 42 d.o.f. If we take a ratio of the flaring and 
the non flaring spectrum, the medium energy band of 6-10 
keV shows the maximum difference. This is in accordance 
with the hardness ratio changes during the flares (Figure [3]). 
The flux measured in the 2-20 keV band during the time 
of the QPO detection is about 1.3xl0 -10 erg cm _2 sec _1 . 
The flux measured in the 2-20 keV band using the spectra 
shown in Figure [6] during the time of flare and non flare time 



Parameter 


Flare 


Nonflare 


QPO 


N H (10 22 cm" 2 ) 


3.96+0.08 


3.03+0.62 


0.283+0.09 


r 


0.98+0.15 


0.79±0.06 


1.55±0.008 


norm 


0.46±0.087 


0.24±0.03 


0.016±0.00045 


E c 


6.29±0.52 


4.56±0.33 


17.07 


Ef 


9.46±0.71 


9.20±0.45 


19.46 


£ G (keV) 


6.57±0.28 


6.56±0.09 


6.42±0.069 


GN 


0.0031±0.0006 


0.0049±0.0009 


0.000086±0.000013 


EQW(keV) 


0.04 


0.12 


0.1 


Xr 


1.37 


1.62 


0.9 



are 7.7x10" 
respectively. 



erg cm sec 



1 and 5.2x10 9 erg cm 2 sec 



4 DISCUSSION 

Detailed timing and spectral studies of 4U 1901+03 using 
this set of RXTE-PCA observations have been reported ear- 
lier. This includes timing and spectral evolution studies over 
the duration of the outburst, pulse frequency evolution (Gal- 
loway et al. 2005), energy and luminosity dependence of the 
pulse profile (Chen et al. 2008) and pulse phase dependence 
of the energy spectrum (Lei et al. 2009). In the previous 
section we have reported some additional temporal studies, 
including the detection of an intensity dependence of the 
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pulse profile at very short time scales; probably for the first 
time in any short period pulsar. 



4.1 Flares 

X-ray Flares of different duration, occurrence and 
flux enhancement have been reported in many accre- 
tion powered puls ars like GRO J174 4-28 (|Finger et all 
Il996h7 SMC X-l (|Moon et all l2003bT). EXO 2030+375 



i Parmar et al.lll9891).LMC X-4 dMoon et al.ll2001bl,l2003a|) . 



4U 1907+ 09 din 't Zand et al.ll 19981; iMukheriee et al.ll200lh . 
Vela X-l jKrevkenbohm et al.ll2008f). SWIFT J1626.6- 5156 
l|Reig et al.ll2008h. Her X-l dMpon fc Eikenberry||2001ah and 
4U 1626-67 ([McClintock et al.lll980h . 

The extremely large flux enhancements seen in GRO 
J1744-28 and the large flares in SMC X-l are probably due 
to accretion disk instability similar to the bursts in the rapid 
burster. The large flares in EXO 2030+375 and LMC X-4 
are similar, and can be caused by Rayleigh- Taylor Instabil- 
ity at the magnetosphere or variable mass accretion from a 
clumpy stellar wind. The flares from the above sources lack 
any associated spectral change. The fast rise and slow decay 
as seen in some of the flares in 4U 1901+03 are similar to 
those in the flares in EXO 2030+375, though the latter were 
of much longer duration. The 1000 s flares from 4U 1626-67 
are caused by carbon burning under the neutron star sur- 
face. The flares in 4U 1901+03 show a spectral hardening, 
while a spectral softening is expected if these flares had been 
of thermonuclear origin. Somewhat smaller amplitude flares 
are seen in 4U 1907+09, Vela X-l and SWIFT J1626. 6-5156 
which are most likely caused by variable wind accretion. 
Overall, the flares in 4U 1901+03 closely resemble the flares 
in SWIFT J1626. 6-5156 in their duration, flux enhancement, 
and spectral variation. It is therefore quite likely that the 
flares in 4U 1901+03 are also caused by variable accretion, 
either due to inhomogeneities in the accretion disk or small 
scale clumping in the stellar wind as is seen in some Super- 
giant Fast X-ray Transients (Sidoli et al. 2010 and references 
therein) . 



Figure 6. The 2.5-25 keV spectra of 4U 1901+03 during X-ray 
flares (top), non-flares (middle) and during QPO detection (bot- 
tom) are shown here along with the respective best fit models con- 
sisting of absorption, power law, high-energy cutoff and gaussian 
components. The lower panel in each figure shows the residual to 
the best fit models. 



4.2 Broadening of the pulse frequency peak and 
intensity dependence of pulse profiles 

A significant broadening of the pulse frequency was dis- 
covered in sources like Cen X-3, 4U 1145-62 and Vela 
X-l using EXOSAT observations. It was proposed that 
the coupling between the periodic and aperiodic variabil- 
ity appears as broadening in the w ings of the spin fre- 
quency peak (jLazzati fe Stella! 1 19971 ). Multiple couplings 
were discovered in Her X-l from RXTE-PCA observations 
(|Moon fe Eikenberrvl I^OOlaT ). Recently, using RXTE-PCA 
observations of several sources, a broaden ing of the pulse 
frequ ency peak has b een shown: Cen X -3 (|Raichur fe Paull 
l2008h. 4U 1626-6 7 l|jain et al.l |2010D and KS 1947+30 
( James et aLll2010h . For most of the outburst period of 4U 
1901+03, we have found a broadening of the pulse frequency 
peak that we attribute to a change in the pulse shape at a 
short time scale of a few tens to a few hundreds of sec- 
onds. Interestingly, the broadening of the pulse frequency 
peak and the pulse shape change are absent in the low state 
when the QPOs appear. 



6 Marykutty James, Biswajit Paul, Jincy Devasia and Kavila Indulekha 



Over the entire out burst period, 4U 1901 +03 showed a 
variety of pulse profiles (|Gallowav et alj|2005l ) . It is interest- 
ing to note that some of the pulse profile changes reported 
in Galloway et al. (2005) at different intensity levels are sim- 
ilar to the profiles shown here, though the former were ob- 
tained at different overall intensity level with time difference 
of several days/weeks and the latter are acquired within a 
timespan of less than one hour. This indicates that the pul- 
sar beaming pattern changes from pulse to pulse, depending 
on the instantaneous mass accretion rate. 



4.3 Quasi Periodic Oscillations 

Quasi Periodic Oscillations from Accretion Powered Pul- 
sars have thus far been reported from 17 sources including 
the recent discovery from the transient pulsar KS 1947+300 
(James et al. 2010). The first discovery of a QPO in an 
accretion powere d pulsar was made i n 1989, in the source 
EXO 2030+375 (|Angelini et alJll989h . All types of accre- 
tion powered pulsars, i.e persistent and transient LMXB & 
HMXB pulsars have shown low frequency QPO features. 4U 
1901+03 is the tenth high magnetic field transient pulsar in 
which QPOs have been detected. 

In accretion powered pulsars, QPOs seem to occur more 
in transient sources compared to the persistent ones. In most 
of these cases, the QPO shows a transient nature. In KS 
1947+300 and 4U 1901+03, the QPOs are detected dur- 
ing the end of the outburst only while a persistent QPO 
feature was seen i n the transient pulsar XTE J1858+034 
ijPaul fc Radll99"^ l. Persi stent pulsars like Cen X-3 do not 
show QPOs all the time llRaichur fc Paulll2008r ) except for 
4U 1626+67 (|Kaur et alj|2008h in which QPOs are detected 
in most of the observations. 

There are many models that have been proposed to 
explain the Quasi Periodic Oscillations in X-ray binaries. 
Beat Frequency Model is generally considered as the most 
plausible model for QPOs in X-ray Pulsars. In this model, 
blobs of matter entrained in the neutron star magnetic field 
which orbit with the Keplerian frequency of the inner edge 
of the accretion disk are accreted at a rate that is modu - 
lated by the rotating magnetic field (|Alpar fc Shahaml 19851 ). 
It produces an aperiodic variability at the beat frequency 
between the pulsar spin frequency and the Keplerian fre- 
quency; vqpo — vu — v s . Another possible model is the 
Keplerian Frequency Model in which the accretion disk con- 
tains structures that persists for a few cycles around the 
neutron star a nd modulate the X-ray f lux by obscuration. 
I.e vqpo = v% ijvan der Klis et al"1ll987f ). 

In the present source and in some other sources like 
Cen X-3, 4U 0115+63, 4U 1626-67, Her X-l, LMC X-4, V 
0332+52, SMC X-l and KS 1947+300, the pulsar frequency 
is higher than the QPO frequency hence the Keplerian fre- 
quency model is inapplicable. It is expected that centrifugal 
forces inhibit mass accretion if the Keplerian frequency at 
the inner disk radius is less than the spin frequency of the 
pulsar. The X-ray flux is proportional to the mass accre- 
tion rate, and the Keplerian frequency at the inn er edge 
of th e accretion disk is related to the X-ray flux l| Finger] 
Il998f l. So in both KFM and BFM we expect a positive cor- 
relation between QPO frequency and the luminosity which 
was however, n ot seen in two of the sources V 0332+5 3 
l|Qu et all [20051 ) and GRO J1744-28 (|Zhang et al.lll996l ). 



So neither of these models are applicable in these sources. 
But in sour ces like A0535+262, E XO 2030+375 and XTE 
J1858+034 l|Mukheriee et al.ll2006l ) the QPO frequency de- 
pends on the X-ray intensity. 

In the 2.7 s pulsar 4U 1901+03, the QPO is detected 
at ~0.135 Hz during decline phase of the 2003 outburst. 
The same kind of QPO feature was also discovered in KS 
1947+300, during its 2001 outburst. KS 1947+300 source 
has the smallest QPO fre quency among all the transient 
pulsars (|james et al]l20ld ). In 4U 1901+03, the QPO fre- 
quency is less than the spin frequency of the pulsar, so it 
can be explained by the Beat Frequency Model. Radius of 
the QPO production region can be estimated as 
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GM 
4tt 2 ^ 



1/3 



(1) 



For the Beat Frequency Model, the Keplerian frequency 
is Uk ~ 0.5 Hz. Assuming a neutron star mass of 1.4 Mq, 
the inner radius of the accretion disk is obtained as 2.6 x 10 3 
km. 

The 2-20 keV X-ray flux during the time of QPO de- 
tection is 1.3xlO~ 10 erg cm -2 sec _1 . The total X-ray lumi- 
nosity for the source is calculated to be about 1.6 x 10 36 erg 
s _1 at 10 kpc. The magnetospheric radius of a neutron star 
can be expr essed in terms of the luminosity and magnetic 
moment as (Fr ank et alJl20oj ). 



Rn 



o , 1 n8 r -2/7 4/7 

3 x 10 L 37 fj, 30 cm 



(2) 



where L37 is the X-ray luminosity in units of 10 erg 
s _1 and /i3o is the magnetic moment in units of 10 30 cm 3 
Gauss. Equating R m with R qpo we estimate the magnetic 
field as 0.31 xlO 12 G. This is consistent with the value in- 
ferred by Galloway et al. (2005) by estimating magnetic field 
strength from spin equilibrium (0.3xl0 12 G). This is also 
within the bounds for the magnetic field derived from the 
upper and lower flux levels within which the source behaved 
as an X-ray pulsar. In the Beat Frequency Model, the ex- 
pected QPO frequency range, for source luminosity between 
5.8 x 10 37 erg s" 1 at the peak of the outburst and 0.16x 10 37 
erg s~ when the QPOs are detected, is 2.05 to 0.135 Hz. 
But only a 0.135 Hz QPO is detected in this source. 



5 ACKNOWLEDGEMENTS 

We thank an anonymous referee for suggestions that helped 
us to improve the manuscript. This research has made use 
of data obtained through the High Energy Astrophysics Sci- 
ence Archive Research Center Online Service, provided by 
the NASA/Goddard Space Flight Center. 



REFERENCES 

Alpar M. A., Shaham J., 1985, Natur, 316, 239 
Angelini L., Stella L., Parmar A. N., 1989, ApJ, 346, 906 
Chen W., Qu J.-L., Zhang S., Zhang F., Zhang G.-B., 2008, 

ChA&A, 32, 241 
Finger, M. H., Koh, D. T., Nelson, R. W., Prince, T. A., 

Vaughan, B. A., Wilson, R. B., 1996, ApJ, 381, 291 



Flares, pulse peak broadening and QPOs in J^U 1901+03 7 



Finger M. H., 1998, Adv. Space Res., 22, 1007 
Forman W., Tananbaum H., Jones C, 1976, ApJ, 206, L29 
Frank J., King A., Raine D. J., 2002, apa..book, 
Galloway D. K., Wang Z., Morgan E. H., 2005, ApJ, 635, 
1217 

in 't Zand, J. J. M., Baykal, A., Strohmayer, T. E. 1998, 

Astroph. J. , 496, 386 
Jain, C, Paul, B., Dutta, A. 2010, MNRAS, 403, 920 
James M., Paul B., Devasia J., Indulekha K., 2010, Ac- 
cepted in MNRAS l|arXiv:1005.116ip 
Kaur, R., Paul, B., Kumar, B., & Sagar, R. 2008, Astroph. 

J. , 676, 1184 
Kreykenbohm, I. et al. 2008, A&A, 492, 511 
Lazzati D., Stella L., 1997, Astroph. J. , 476, 267 
Lei Y.-J., Chen W., Qu J.-L., Song L.-M., Zhang S., Lu Y., 

Zhang H.-T., Li T.-P, 2009, Astroph. J. , 707, 1016 
McClintock, J. E., Canizares, C. R., Li, F. K., & Grindlay, 

J. E. 1980, Astroph. J. , 235, L81 
Molkov S., Lutovinov A., Grebenev S., 2003, A&A, 411, 
L357 

Moon D.-S., Eikenberry S. S., 2001a, Astroph. J. , 552, 
L135 

Moon D.-S., Eikenberry S. S., 2001b, Astroph. J. , 549, 
L225 

Moon D.-S., Eikenberry S. S., Wasserman I. M., 2003a, 

Astroph. J. , 586, 1280 
Moon D.-S., Eikenberry S. S., Wasserman I. M., 2003b, 

Astroph. J. , 582, L91 
Mukherjee K., Agrawal, P. C, Paul, B., Rao, A. R., Yadav, 

J. S., Seetha, S., Kasturirangan, K. 2001, Astroph. J. , 

548, 368 

Mukherjee U., Bapna S., Raichur H., Paul B., Jaaffrey 

S. N. A., 2006, JApA, 27, 25 
Parmar A. N., White N. E., Stella L., Izzo C, Ferri P., 

1989, Astroph. J. , 338, 359 
Paul B., Rao A. R., 1998, A&A, 337, 815 
Priedhorsky W. C, Terrell J., 1984, Astroph. J. , 280, 661 
Qu J. L., Zhang S., Song L. M., Falanga M., 2005, ApJ, 

629, L33 

Raichur H., Paul B., 2008, Astroph. J. , 685, 1109 
Reig P., Martmez-Nunez S., Reglero V., 2006, A&A, 449, 
703 

Reig, P., Belloni, T., Israel, G. L., Campana, S., Gehrels, 

N., Homan, J., 2008, A&A, 485, 797 
Sidoli, L. Esposito, P. Ducci, L., 2010, Mon. Not. Roy. Ast. 

Soc. (in press), astro-ph: 1007. 1091 
van der Klis M., Stella L., White N., Jansen F., Parmar 

A. N., 1987, Astroph. J. , 316, 411 
Zhang W., Morgan E. H., Jahoda K., Swank J. H., 

Strohmayer T. E., Jernigan G., Klein R. I., 1996, ApJ, 

469, L29 



